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In a companion paper (see Resource Allocation for Downlink Cellular OFDMA Systems: Part I — 
Optimal Allocation), we characterized the optimal resource allocation in terms of power control and 
^ ' subcarrier assignment, for a downlink sectorized OFDMA system impaired by multicell interference. In 

our model, the network is assumed to be one dimensional (linear) for the sake of analysis. We also 
^f) . assume that a certain part of the available bandwidth is likely to be reused by different base stations 

■ while that the other part of the bandwidth is shared in an orthogonal way between these base stations. 
The optimal resource allocation characterized in Part I is obtained by minimizing the total power spent 
by the network under the constraint that all users' rate requirements are satisfied. It is worth noting that 
when optimal resource allocation is used, any user receives data either in the reused bandwidth or in the 
protected bandwidth, but not in both (except for at most one pivot-user in each cell). We also proposed 
an algorithm that determines the optimal values of users' resource allocation parameters. 
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, As a matter of fact, the optimal allocation algorithm proposed in Part I requires a large number of 

^ . operations. In the present paper, we propose a distributed practical resource allocation algorithm with low 

complexity. We study the asymptotic behavior of both this simplified resource allocation algorithm and 
the optimal resource allocation algorithm of Part I as the number of users in each cell tends to infinity. 
Our analysis allows to prove that the proposed simplified algorithm is asymptotically optimal i.e., it 
achieves the same asymptotic transmit power as the optimal algorithm as the number of users in each 
cell tends to infinity. As a byproduct of our analysis, we characterize the optimal value of the frequency 
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reuse factor. Simulations sustain our claims and show that substantial performance improvements are 
obtained when the optimal value of the frequency reuse factor is used. 

Index Terms 

OFDMA, Multicell Resource Allocation, Distributed Resource Allocation, Asymptotic Analysis. 

I. Introduction 

In a companion paper HI, we introduced the problem of joint power control and subcarrier assignment 
in the downUnk of a one-dimensional sectorized two-cell OFDMA system. Resource allocation parameters 
have been characterized in such a way that i) the total transmit power of the network is minimum and 
ii) all users' rate requirements are satisfied. Similarly to |2|, we investigate the case where the channel 
state information at the Base Station (BS) side is limited to some channel statistics. However, contrary 
to JH, our model assumes that the available bandwidth is divided into two bands: the first one is reused 
by different base stations (and is thus subject to multicell interference) while the second one is shared in 
an orthogonal way between the adjacent base stations (and is thus protected from multicell interference). 
The number of subcarriers in each band is directly related to the frequency reuse factor. We also assume 
that each user is likely to modulate subcarriers in each of these two bands and thus we do not assume a 
priori a geographical separation of users modulating in the two different bands. The solution to the above 
resource allocation problem is given in the first part of this work. This solution turns out to be "binary": 
except for at most one pivot-user, users in each cell must be divided into two groups, the nearest users 
modulating subcarriers only in the reused band and the farthest users modulating subcarriers only in the 
protected band. An algorithm that determines the optimal values of users' resource allocation parameters 
is also proposed in the first pait. 

It is worth noting that this optimal allocation algorithm is still computationally demanding, especially 
when the number of users in each cell is large. One of the computationally costliest operations involved 
in the optimal allocation is the determination of the pivot-user in each cell. In the present paper, we 
propose a distributed simplified resource allocation algorithm with low computational complexity, and 
we discuss its performance as compared to the optimal resource allocation algorithm of Part I. This 
simplified algorithm assumes a pivot-distance that is fixed in advance prior to the resource allocation 
process. Of course, this predefined pivot-distance should be relevantly chosen. For that sake, we show 
that when the fixed pivot-distance of the simplified algorithm is chosen according to a certain asymptotic 
analysis of the optimal allocation scheme, the performance of the simplified algorithm is close to the 
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optimal one, provided that the number of users in the network is large enough. Therefore, following 
the approach of we propose to characterize the limit of the total transmit power which results from 
the optimal resource allocation policy as the number of users in each cell tends to infinity. Several 
existing works on resource allocation resorted to this kind of asymptotic analysis, principally in order to 
get tractable formulations of the optimization problem that can be solved analytically. For example, the 
asymptotic analysis was used in [31 and |4J in the context of downlink and uplink single cell OFDMA 
systems respectively, as well as in Q in the context of Code Division Multiple Access (CDMA) systems 
with fading channels. Another application of the asymptotic analysis can be found in 0. The authors 
of the cited work addressed the optimization of the sum rate performance in a multicell network. In this 
context, the authors proposed a decentralized algorithm that maximizes an upper-bound on the network 
sum rate. Interestingly, this upper-bound is proved to be tight in the asymptotic regime when the number 
of users per cell is allowed to grow to infinity. However, the proposed algorithm does not guaranty 
fairness among the different users. 

In this paper, we use the asymptotic analysis in order to obtain a compact form of the (asymptotic) 
power transmitted by the network for the optimal resource allocation algorithm, and we use this result 
to propose relevant values of the fixed pivot-distance associated with the simplified allocation algorithm. 
We prove in particular that when this fixed pivot-distance is chosen equal to the asymptotic optimal 
pivot-distance, then the power transmitted when using the proposed simplified resource allocation is 
asymptotically equivalent to the minimum power associated with the optimal algorithm. This limiting 
expression no longer depends on the particular network configuration, but on an asymptotic, or "aver- 
age", state of the network. More precisely, the asymptotic transmit power depends on the average rate 
requirement and on the density of users in each cell. It also depends on the value a of the frequency 
reuse factor. As a byproduct of our asymptotic analysis, we are therefore able to determine an optimal 
value of the latter reuse factor. This optimal value is defined as the value of a which minimizes the 
asymptotic power. 

The rest of this paper is organized as follows. In Section |ll] we recall the system model as well as 
the joint resource allocation problem. In Section |llll we propose a novel suboptimal distributed resource 
allocation algorithm. Section |IV] is devoted to the asymptotic analysis of the performance of this simplified 
allocation algorithm as well as the performance of the optimal resource allocation scheme of Part I when 
the number of users tends to infinity. Theorem [T] characterizes the asymptotic behavior of the optimal 
joint allocation scheme. The results of this theorem are used in Subsection IIV-DI in order to determine 
relevant values of the fixed pivot-distances associated with the simplified allocation algorithm. Provided 
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that these relevant values are used, Proposition |2] states that the simplified algorithm is asymptotically 
optimal. Section |Vl] addresses the selection of the best frequency reuse factor. Finally, Section IVIII is 
devoted to the numerical illustrations of our results. 



A. System Model 

We consider a sectorized downlink OFDMA cellular network. We focus on two neighboring one- 
dimensional (linear) cells, say Cell A and Cell B, as illustrated by Figure [U Denote by D the radius of 



Figure 1. Two-Cell System model 

each cell. We denote by the number of users of Cell A and by the number of users of Cell B. The 
total number of available subcarriers in the system is denoted by N. For a given user k G 1,2, ... , K'^ 
in Cell c (c G {A, B}), we denote by Xk the distance that separates him/her from BS c, and by ]\ffc the 
set of indices corresponding to the subcarriers modulated by k. J^k is a subset of {0, 1, . . . , — 1}. The 
signal received by user k at the nth subcarrier (n G Dsf^) and at the mth OFDM block is given by 



where Sk{n,'m) represents the data symbol transmitted by BS c. Process Wk{n,m) is an additive noise 
which encompasses the thermal noise and the possible multicell interference. Coefficient Hk{n,m) is 
the frequency response of the channel at the subcarrier n and the OFDM block m. Random variables 
Hf^{n,m) are assumed Rayleigh distributed with variance = E[\Hk{n,m)\'^]. Channel coefficients 



II. System Model and Previous Results 



Frequency 




yk{n, m) = Hk{n, m)sk{n, m) + Wk{n, m). 



(1) 
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are supposed to be perfectly known at the receiver side, and unknown at the BS side. We assume that 
Pk vanishes with the distance Xk based on a given path loss model. The set of available subcarriers is 
partitioned into three subsets: J containing the reused subcarriers shared by the two cells; and 
containing the protected subcarriers only used by users in Cell A and B respectively. The reuse factor 
a is defined as the ratio between the number of reused subcarriers and the total number of subcarriers: 

card(J) 
" = ^^ 

so that 3 contains aN subcarriers. If user k modulates a subcarrier n G J, the additive noise contains both 
thermal noise of variance and interference. Therefore, the variance a\ of this noise-plus-interference 



process depends on k and coincides with cr? = E 



\Hk{n,m) 



Qf + fj^, where Hk{n,m) represents 



the channel between BS B and user k of Cell A at frequency n and OFDM block m, and where Qf = 
^ki^ki average power transmitted by BS B in the interference bandwidth J. The remaining 

{\ — Q.)N subcarriers are shared by the two cells. Cell A and i? , in an orthogonal way. If user k modulates 
such a subcarrier n G !Pc, the additive noise Wk{n,m) contains only thermal noise. In other words, 
subcarrier n does not suffer from multicell interference. Then we simply write E,[\wk{n, = o"^. The 
resource allocation parameters for user k are: P| ^ the power transmitted on each of the subcarriers of 
the non protected band 3 allocated to him, 7^ ^ his share of 3, Pj: ^ the power transmitted on each of the 
subcarriers of the protected band allocated to him and 7^ ^ his share of Tc. In other words, 

7^,1 = card(J n Jik)lN 7^^2 = n Ji^)lN . 

As a consequence, J2^=i lki~'-^ Ylik=i 7fc 2 ~ ^'^^ e-ach cell c. Moreover, let g^^i (resp. gk^2) 
be the channel Gain to Noise Ratio (GNR) in band 3 (resp. IPc), namely g^ i = Pk/c^'i (resp. gk^2 = Pfc/c^). 
"Setting a resource allocation for cell c" means setting a value for parameters {7^ ii 7fc 2' ^fc i' Pk2\k=i...K''- 

B. Joint Resource Allocation for Cells A and B 

Assume that each user k has a rate requirement of Rk nats/s/Hz. In the first Part of this work (H, our 
aim was to jointly optimize the resource allocation for the two cells which i) allows to satisfy all target 
rates Rk of all users, and ii) minimizes the power used by the two base stations in order to achieve these 
rates. For each cell c € {A, B}, denote by c the adjacent cell {A = B and B = A). The ergodic capacity 
associated with a user k in Cell c is given by 

Ck = 7mIE [log (1 + gk,imPk,iZ)] + 7fc,2lE [log (l + 5fe,2n^,2^)] , (2) 
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where Z is a standard exponentially distributed random variable, and where coefficient gk,i{Q'i) is given 
by 

OkAQl) = , (3) 



E \Hk{n,m)\-^ Qi + 

where i7fc(n, m) represents the channel between BS c and user k of Cell c at frequency n and OFDM 
block m. Coefficient gk,i{Qi) represents the signal to interference plus noise ratio in the interference 
band J. We assume that users are numbered from the nearest to the BS to the farthest. As in [l], 
the following problem will be refeiTcd to as the joint resource allocation problem for Cells A and B: 

Minimize the total power spent by both base stations = ^^(7A;,i-Pfc,i +7fc,2^fc,2) with respect 

c=A,Bk=l 

{7fc 1'7a: 2'-^! I'-^fc 2} c=A,B Under the following constraint that all users' rate requirements are 

' k=l..'.K'= 

satisfied i.e., for each user k in any cell c, Rk < Ck- The solution to this problem has been determined in 
the first part of this work |fT1|. As a noticeable point, the results of 11] indicate the existence in each cell 
of a pivot-user that separates two groups of users: the "protected" users and the "non protected" users. 
The following proposition states this binary property of the solution. 

Proposition 1 ([1|). Any global solution to the joint resource allocation problem is "binary" i.e., there 
exists a user L'^ in each Cell c such that 7/^ 2 = far closest users k < L^, and 7/^1 = far farthest 
users k > L'^. 

In the sequel, we denote by d^'^^^ the position of the pivot-user in Cell c i.e., d'^'^^^ = xl"- A 
resource allocation algorithm is also proposed in [H. This algorithm turns out to have a high computational 
complexity and the determination of the optimal value of the pivot-distance d^'^^-* turns out to be one 
of the costliest operations involved in this algorithm. This is why we propose in the foUwing section of 
the present paper a suboptimal simplified allocation algorithm that assumes a predefined pivot-distance. 

III. Practical Resource Allocation Algorithm 

A. Motivations and Main idea 

Proposition [T] provides the general form of the optimal resource allocation, showing in particular the 
existence of pivot-users L"^, in both Cells A, B, separating the users who modulate in band J from the 
users who modulate in bands 7^ and 7^ . As a matter of fact, the determination of pivot-users L^, is 
one of the costliest operations of this optimal allocation (see [l] for a detailed computational complexity 
analysis). Thus, it would be convenient to propose an allocation procedure for which the pivot-position 
would be fixed in advance to a constant rather than systematically computed/optimized. We propose a 
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simplified resource allocation algorithm based on this idea. Furthermore, we prove that when the value 
of the fixed pivot-distances is relevantly chosen, the proposed algorithm is asymptotically optimal as the 
number of users increases. In other words, the total power spent by the network for large K when using 
our suboptimal algorithm does not exceed the minimum power that would have been spent by using the 
optimal resource allocation. The proposed algorithm is based on the following idea. 

Recall the definition of (i^'(^) and d^'^^^^ as the respective position of the optimal pivot-users and 
defined by Proposition [2 As the optimal pivot-positions d'^'^^^ and d^'^^^ are difficult to compute 
explicitly and depend on the particular rates and users' positions, we propose to replace d-^'^^^ and d^'(^) 
with predefined values (i^ibopt (i^bopt fixed before the resource allocation process. In our suboptimal 
algorithm, all users in Cell c whose distance to the BS is less than dsubopt modulate in the interference 
band 3. Users farther than (isubopt modulate in the protected band J"^. Of course, we still need to determine 
the pivot-distances d^^^^^ and (i^bopt- ^ procedure that permits the relevant selection of d^bopt? <^^bopt 
is given in Section IIV-CI 



B. Detailed Description 

Assume that the values of fi^bopt ^^'^ "^Sbopt ^^^^ been fixed beforehand prior to the resource allocation 
process. For each Cell c, define by %'j the subset of {1, . . . K'^} corresponding to the users whose distance 
to BS c is less than djubopt- Define by Xp the set of users whose distance to BS c is larger than dsubopt- 

1) Resource allocation for protected users: Focus for instance on Cell A. For each k G 3Cp, we 
arbitrarily set = P^^ = i.e., user k is forced to modulate in the protected band only. For 
such users, the remaining resource allocation parameters 7^2) Pk2 obtained by solving the following 
classical single cell problem w.r.t. {jk2^-^k'2)kex^' 

"Minimize the transmitted power Ylk^x^ 1k2-^k2 under rate constraint Rk < Ck for each k G Xp". 

The above problem is a simple particular case of the single cell problem addressed in fV\. Define the 
functions f{x) = — ^ and C{x) = E[log(l + f^^{x)Z)] on M_|_. The solution is given by 

Pk^,2 = 9kirH9k,2P2) 

A 



E 



log (1+5^,2^^2^ 
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where parameter [32 is obtained by writing that constraint = holds or equivalently, (32 is 

the unique solution to: 

\ - Rk _ 1- a 

,%c{a,,2h)~ 2 

We proceed similarly for Cell B. 

2) Resource allocation for interfering users: We now focus on users k € for each cell c = A,B. For 
such users, we arbitrarily set 7^ 2 = -ffc 2 ~ *-* users in Xj are forced to modulate in the interference 
band J only, for each cell c. The remaining resource allocation parameters ^ki->-^ki->lki^ ^ki obtained 
by solving the following simplified multicell problem. 

Problem 1. [Multicell] Minimize ^ ^ 7fc,i-Pfc,i {lki^^ki^^ki^^ki)k under the following 

c=A,Bk(^OC'=j 

constraints for each cell c G {A,B}: 

CI : Vc, Vfc G Xj, Rk <Ck C2 : Vc, ^ Yk,i = « C3 : -fk,i > • 

kf^OCj 

Clearly, the above Problem can be interpreted as a particular case of the initial resource allocation 
(Problem 2 in [ID addressed in Section Hl-B I of the present paper. The main difference is that the initial 
multicell problem jointly involves the resource allocation parameters in three bands J, "J"^ and whereas 
the present problem only optimizes the resource allocation parameters corresponding to band J, while 
arbitrarily setting the others to zero. Therefore, the results of Part I ID, Theorem 2 of ID in particular, 
can directly be used to determine the global solution to Problem [D 

Remark 1 (Feasibility). Recall that the initial joint resource allocation Problem (Problem 2 in [T]) 
described in Section \II-B\ in the present paper was always feasible. Intuitively, this was due to the fact 
that any user was likely to modulate in the protected band if needed, so that any rate requirement Rk 
was likely to be satisfied by simply increasing the power in the protected band. In the present case, the 
protected band is by definition forbidden to users in Xp Theoretically speaking. Problem \J} might not be 
feasible due to multicell interference. Fortunately, we will see this case does not happen, at least for a 
sufficiently large number of users, if the values of the pivot-distances d^j^^p, and df^/^^^^ are well chosen. 
This point will be discussed in more detail in Section [3 

Define = J^keX" Tfc i^k 1 average power transmitted by BS c in the interference bandwidth 3. 
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By straightforward application of Theorem |2j we obtain that for each Cell c and for each user k € JCj, 

Pk,i=9;liQl)r\9kAQlWi) (4) 
111 = XT , (5) 



where for each c = A,B and for a fixed value of , parameters Ql) are the unique solution to the 
following system of equations: 

y = a (6) 

Ui = > Kk z — = . (/) 

Note that the first equation is nothing else that the constraint C2: ^^7^ ^ = a. The second equation 
is nothing else than the definition = Ylki^oc i-^k v prove that the system of four 

equations ©-([T]) for c = A,B admits a unique solution (3i , Qf , (3^ , Qf and we provide a simple 
algorithm allowing to determine this solution. 

Focus on a given Cell c and consider any fixed value Q^. Denote by /^(Qf ) the rhs of equation ^ 
where /3J is defined as the unique solution to ®. Since ([7]) should be satisifed for both c = A and c = B, 
the following two equations hold 

Qt = i^iQf), Qf = i''iQt)- 

The couple {Qf,Qi) is therefore clearly a fixed point of the vector-valued function I{Qf',Qf) = 
{i^{Qf),P{Qt)). 

{QtQf) = iiQt,Qf)- (8) 

As a matter of fact, it can be shown that such a fixed point of I is unique. This claim can be proved 
using the approach previously proposed by [12]. 

Lemma 1. Function I is such that the following properties hold. 

1) Positivity: i{Q^,Q^) > 0. 

2) Monotonicity: If > Q'^'.Q^ > Q^', then i{Q^,Q^) > i{Q^',Q^'). 

3) Scalability: for all t > 1, ti{Q^,Q^) > l{tQ^,tQ^). 

The proof of Lemma [T] uses arguments which are very similar to the proof of Theorem 1 in [11]. It is 
thus omitted from this paper and provided in |[T3l . Function I is then a standard interference function, 
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using the terminology of |[T2l . Therefore, as stated in |[T2l . such a function I admits at most one fixed 
point. On the other hand, the existence of a fixed point is ensured by the feasibiUty of Problem [T] and by 
the fact that dS) holds for any global solution. In other words, if Problem \T\ is feasible, then function I 
does admit a fixed point and this fixed point is unique. Putting all pieces together, there exists a unique 
solution to ([8]), which can be obtained thanks to a simple fixed point algorithm. In practice, resource 
allocation in band 3 can be achieved by the following procedure. 
Ping-pong algorithm for interfering users 

1) Initialization: Qf = 0. 

2) Cell A: Given the current value of the power Qf transmitted by base station B in the interference 
bandwidth, compute pf, Qf- as the unique solution to ©-(IT]) with c = A. 

3) Cell B: Given the current value of Qf, compute /3f , Qf by ©-d?]). 

4) Go back to step 2 until convergence. 

5) Define resource allocation parameters by ©-([S]). 
Comments 

1) Convergence of the ping-pong algorithm. We stated earlier that Problem [T] is either feasible or 
infeasible, depending on the value of (d^i^^pp df^^^^p^). If the latter problem is feasible, then function 
i will heve a unique fixed point due to Lemma [T] and the ping-pong algorithm will converge to 
this fixed point. If Problem \T\ is infeasible, then function I will have no fixed points and the the 
ping-pong algorithm will diverge. One of the main purposes of Section IIV-CI is to provide relevant 
values of (d^bopf '^Sbopt) ^^^^ that convergence of the ping-pong algorithm holds for sufficiently 
large number K of users. 

2) Note that the only information needed by Base Station c about Cell c is the current value of 
the power Ql transmitted by Base Station c in the interference band J. This value can i) either 
be measured by Base Station c at each iteration of the ping-pong algorithm, or ii) it can be 
communicated to it by Base Station c over a dedicated link. In the first case, no message passing 
is required, and in the second case only few information is exchanged between the base stations. 
The ping-pong algorithm can thus be implemented in a distributed fashion. 

C. Complexity Analysis 

We showed earlier that allocation for protected users can be reduced to the determination in each cell 
of the value of /Jg, which is the unique solution to the equation ^i-^x-'^ ~ argued 

in IT] that solving this kind of equations requires a computational complexity proportional to the number 
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of terms in the Ihs of the equation, which is itself of order 0{K). Using similar arguments, we can 
show that each iteration of the ping-pong algorithm for non protected users can be performed with a 
complexity of order 0{K). Let J designate the number of iterations needed till convergence. The overall 
computational complexity of the ping-pong algorithm, and hence of the simplified resource allocation 
scheme as well, is thus of the order of 0{JK). Our simulations showed that the ping-pong algorithm 
converges relatively quickly in most of the cases. Indeed, no more than J = 15 iterations were needed in 
almost all the simulations settings to reach convergence within a very reasonable accuracy. The complexity 
of the simplified algorithm is to be compared with the computational complexity of the optimal algorithm 
which was shown in lU to be of the order of 0{MK \og2 K), where M is the number of points inside 
a certain 2D search grid. 

IV. Asymptotic Optimality of the Simplified Resource Allocation Scheme 

The aim of this section is to evaluate the performance of the proposed simplified algorithm. The 
relevant performance metric in the context of this paper is the total power that must be transmitted by 
the base stations. Since the simpUfied algorithm assumes predefined pivot-distances (d^bop,, c^^bopt) fixed 
prior to the resource allocation process, the performance of the proposed algorithm depends on the choice 
of these fixed pivot-distances. One must therefore determine what relevant value should be selected for 
(^sJibopf ^Sbopt)- ^ possible method is addressed in this section and consists in studying the case where 
the number of users tends to infinity. 

A. Main Tools: Asymptotic analysis 

We study first the performance of the optimal allocation algorithm proposed in Part I [1] when the 
number of users in each cell tends to infinity. From the results of this asymptotic study, we conclude the 
asymptotic behaviour of the optimal pivot-distances [d^'^^\ d^'^^^). It turns out that when the number 
K of users increases, the optimal pivot-distances as well as the total transmitted power no longer depend 
on the particular cell configuration, but on an asymptotic state of the network, such as the average rate 
requirement and the density of users in each cell. Thanks to this result, we can now choose the fixed 
pivot-distances associated with the simplified algorithm to be equal to the asymptotic pivot-distances. 
In this case, one can show that the performance gap between the simplified and the optimal allocation 
schemes vanishes for high numbers of users. We introduce now the mathematical assumptions and tools 
that we use for defining the asymptotic regime. 
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1 ) Notations and Basic Assumptions: In the sequel, we denote by B the total bandwidth of the system 
in Hz. We consider the asymptotic regime where the number of users in each cell tends to infinity. We 
denote by = BR]^ the data rate requirement of user k in nats/s, and we recall that R]^ is the data rate 
requirement of user k in nats/s/Hz. Notice that the total rate 'Yl,k=i '''k which should be delivered by BS c 
tends to infinity as well. Thus, we need to let the bandwidth B grow to infinity in order to satisfy the 
growing data rate requirement. Recalling that K = + denotes the total number of users in both 
cells, the asymptotic regime will be characterized by ^ oo, i? ^ cxd and K/B — > t where t is a positive 
real number. We assume on the other hand that K'^jK (c G {^,5}) tends to some positive constant as 
K tends to infinity. Without restrictions, this constant is assumed in the sequel to be equal to 1/2 i.e., 
the number of users becomes equivalent in each cell. In order to simplify the proofs of our results, we 
assume without restriction that for each k, the rate requirement is upper-bounded by a certain constant 
^max> < ^max» whcrc Tmax Can be choscn as large as needed, and that users of each cell are located 
in the interval [e, D\ where e > can be chosen as small as needed. Recall that Xk denotes the position 
of each user k i.e., the distance between the user and the BS. The variance of the channel gain of user 
k will be written as pk = p{xk) where p{x) models the path loss. Typically, function p{x) has the form 
p{x) = Xx~^ where A is a certain gain and where s is the path-loss coefficient, s > 2. In the sequel, we 
denote by g2{x) = the received gain to noise ratio in the protected bandwidth, for a user at position 
X. This way, g2{xk) = gk,2- Similarly, we define for each user k in cell A, gi{xk,Qf) = gk,i{Qi)- 
More generally, gi{x,Q) denotes the gain-to-interference-plus-noise ratio in the interference bandwidth 
at position x when the interfering cell is transmitting with power Q in band J. Functions gi{x,.) and 
g2{x) are assumed to be continuous functions of x. It is worth noting that for each x, g2{x) = gi{x,Q). 
Finally, recall that coefficient 7^ ^ (resp. 7^ 2) is defined as the ratio between the part of the interference 
bandwidth J (resp. protected bandwidth !Pc) and the total bandwidth. Thus, 7^ ^ and 7^ 2 tend to zero as 
the total bandwidth B tends to infinity for each k. 



2) Statistical Tools and Main Ideas of the Asymptotic Study: Theorem 2 of Part I HI reduces the 
determination of the whole set of resource allocation parameters in both cells to the determination of 
ten unknown parameters {Qi, Pf, L^,C'^}c=A,B,i=i,2- Parameter in particular represents the power 
transmitted by Cell c in the non protected band 3. Consider now one of the two Cells c € {A, B}, and 
denote by c the second (adjacent) cell. In the sequel, we use the notation Q-^^ (resp. Q2 ) instead of 
Q1 (resp. Q2) to designate the power transmitted by BS c in the non protected band 3 (resp. the protected 
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band Tc) when the optimal solution characterized by Proposition [T] is used. 

Q"/""^ = Etm^m (9) 

k=l 

Qf""^ = Y.^hPi2- (10) 

fc=L= 

The new notation Q^i^^\ Q2^^^ is used to indicate the dependency of the results on the number of users 
K. For the same reason, parameters L"^, /3f , /Jj, will be denoted in the sequel by L'^'^^\ f3i^^\ ff^^^^ , 
^'^^K) respectively. Our goal now is to characterize the behavior of the resource allocation strategy as 
K,B ^ CO and, in particular, the behavior of powers Qi^^\ Q2^^^- By straightforward application of 
Theorem 2 of Part I, Ql'^^^ = Ylk=i Ik i^ki '^^'^ written as 

Q?'^''^= E RkHxk,fi'i^''\Qf^^''\e^^''^) + wi,^,^,, (11) 

where W^^^^^) i = 7£c,(k) i-ffccif) i denotes the power transmitted to the pivot-user L'^'(^) in the inter- 
ference band 3, and where function 3" is defined by 

for each x, /5, Q. The first term in the rhs of (fTTl) represents the total power allocated to all users k < L^'^^\ 
It is quite intuitive that the power allocated to one user W^^^^k) i is negligible when compared to the 
power allocated to all users k < L'^'^^\ Indeed, it will be shown in Appendix |A] that the first term of ([TT]) 
is bounded as K ^ oo wheras VF£e,(K) i tends to zero. In the sequel, we use notation W^^^^^) i = ok{^), 
where ox(l) stands for any term which converges to zero as K ^ oo. In order to study the limit of 
this expression as K tends to infinity, we introduce for each one of the two cells the following measure 
defined on the Borel sets of M_|_ x M_|_ as follows 

1 

U^'i^){I,J) = —Y^6r,,.,{I,J) (13) 
k=l 

where / and J are any intervals of and where dr^,Xk is the Dirac measure at point (rfc,Xfc). In order 
to have more insights on the meaning of this tool, it is useful to remark that i/^'^^^I, J) is equal to 

c {K)^j j\ number of users located in J and requiring a rate (in nats/s) in interval I 
' total number of users 

Thus, measure zy'='(^) can be interpreted as the distribution of the set of couples {rk,Xk) of Cell c. The 

introduction of the above measure simplifies considerably the asymptotic study of the transmit power. 
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Indeed, replacing Rk (in nats/s/Hz) by ^^^^^^^ in equation (ITTI ). we obtain 

where integration is considered with respect to the set A\'^^^ = [0,rmax] x [e, d^'^^^], where d'^'(^) = 
xia,{K) is the position of pivot-user L'^'^^^ and where e can be chosen, as stated earlier in this section, as 
small as needed. It is quite intuitive that the asymptotic power limi^'.^oo Qi^^^ can be obtained from (fT4l) 
by replacing '^ = §X'^bytxi and the distribution ly^'^^^ by the asymptotic distribution of 
couples {rk,Xk) as K tends to infinity. The existence and the definition of this asymptotic distribution 
is provided by the following assumption. 

Assumption 1. As K tends to infinity, measure zy'^'(^) converges weakly to a measure v^. 

We refer to fT\ for the materials on the convergence of measures. In order to have some insight 
on the behavior of equation ([141 ) in the asymptotic regime, imagine for the sake of simplicity that 
sequences d^'(^), d^' Qf '^^^ Qf /jf /jf ^^'W, e^'(^) are convergent and that they 
converge respectively to d^,d^ ,Qf, Qf , P^, , . This assumption is of course arbitrary for the 

moment, but it allows to better understand the main ideas of our asymptotic analysis. More rigorous 
considerations on the convergence of these sequences will be discussed later on. Ignoring at first such 
technical issues, it is intuitive from equation (fT4l ) that Q^'^^^ converges to a constant defined by 

Qi = \jj^^ r-3-(x, /3J, Ql e)dy'i.r, x) , (15) 

where = [0,rjnax] x [ei^^]- In other words, we manage to express the limit of the power Ql'^^^ 
transmitted by station c in the interference band as a function of the asymptotic cell configuration. In 
order to further simplify the above expression, it is also realistic to assume that measure z/'^ is the measure 
product of a limit rate distribution times a limit location distribution. Assumption |2] below is motivated 
by the observation that in practice, the rate requirement of a given user is usually not related to the 
position Xk of the user in each cell. 

Assumption 2. Measure v'^ is such that di''^{r,x) = dC,^{r) x d\^{x) where C,^ is the limit distribution 
of rates and is the limit distribution of the users' locations. Here x denotes the product of measures. 

Measures C, and A respectively correspond to the distributions of the rates and the positions of the users 
within one cell. For instance, the value r"^ = ^ /o^"""' dC,^{r) represents the average rate requirement per 
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channel use in Cell c. We furthermore assume that measures A"^ and are absolutely continuous with 
respect to the Lebesgue measure on [e,D]. Using Assumption |2l equation (ITSl) becomes 

Ql = j^^ nx, e) d\\x). (16) 

Of course, a similar result can be obtained for power transmitted by base station c in 



the protected band J"^. To that end, we simply note that function g2{x) satisfies g2{x) = gi{x,0). Using 



similar tools, the expression of Q^2^^ given by (l24l ) converges as K — > cxo toward 



Q\ = f^ r 'S{x,(ilMd\\x). (17) 

Equations (fT6l ) and ([TT] ) respectively provide the limits of Q\^^^ and Q'^^^^^ ^ function of some 
parameters d'^,Pf,(32 and (assumed for the moment to be the limits of d'^'^^\ (3l'^^\ j32^^^ and 
QI'^^^ as long as such limits exist). These unknown parameters still need to be characterized. Therefore, 
we must determine a system of equations which is satisfied by these parameters. This task is done by 
Theorem [T] given below. 

B. Asymptotic Performance of the Optimal Resource Allocation 

Define the following function ^{x, /?, Q, ^) = / ^ for each x, Q, ^. The proof of the following 

result is provided in Appendix [A] 

Theorem 1. Assume that K = + oo in such a way that K/B ^ t > and jK 1/2. 

Assume that the optimal solution for the joint resource allocation problem {Problem 2 in /[7]/j is used for 
each K. The total power spent by the network Q^'^ = X]c=A B Z]fc^i(7fc i-^fc i +7fc 2^fc 2) converges to 
a constant Qt- The limit Qt has the following form: 

qt = E f r '^1'' "^^'("^^ + r 0' 0) ^^'(^)) ' (IS) 

c=A,B ^-^^ ^ 

where for each c = A,B, the following system of equations in variables d'^, P^, (32,S,^ is satisfied: 

g(x, (31, Qle) dX'ix) = a (19) 

£ 9{x, /32^ 0, 0) dX^ix) = (20) 
^il^^ (^i|ff^/??) = 9.idlFig2idlPI) (21) 
f'^ Qle) dX'ix) = Ql . (22) 
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Moreover, for each c = A,B and for any arbitrary fixed value {Qf'jQi), the system of equations 
(|79l)-(|20l)-(|27|)-(|22]) admits at most one solution {d", P^, P^,^''). 

As a consequence, when optimal multicell resource allocation is used, the total power spent by the 
network converges to a constant which can be evaluated through the results of Theorem [T] This result 
allows to evaluate the asymptotic power spent by the network as a function of the reuse factor a, the 
average rate requirement f and the asymptotic distribution of users in each cell A. 

Now that the asymptotic performance of the optimal allocation scheme has been studied, the value of 
the fixed pivot-distances (i^bopf ^Sbopt associated with the simplified allocation algorithm can be relevantly 
chosen to be equal in each Cell c to the asymptotic pivot distance d'^ defined by Theorem [T] 

C. Determination of the fixed pivot-distances df^ibopf '^fuhopt f^^ simplified allocation scheme 

We stated earlier in Section |lll] that the suboptimal algorithm replaces the optimal value (i^'(^) of 
the pivot-distance in each Cell c with a fixed value dsubopt- Intuitively, if (i^ibopt and (i^ibopt are chosen 
such that d"^'^^^ ~ ^^bopt ^'^^ d^'^^^ ~ ^^bopt f^'" large K, the performance of our algorithm shall be 
close to the optimal one as K increases. Therefore, we must determine an asymptotically optimal pair 
of pivot-distances {d^,d^). To that end we propose the following procedure. 

Note first by referring to Theorem [T] that the value of d^, d^ can be easily determined once the relevant 
values of and Qf have been determined. The remaining task is thus the determination of the value 
of {Qt^ Qi )- To that end, we propose to perform an exhaustive search on {Qi^, Qf )- 
i) For each point {Qf, Qf) on a certain 2D search grid, solve the system ([T9l)- (l20l) -(l2TI)-(l22l) introduced 
by Theorem [T] for both a = A,B. Theorem [T] states that this system admits at most one solution for any 
arbitrary fixed value {Qf, Qf)- If the investigated point {Qf, Qf) of the grid is such that the system (fT9l )- 
(l20l)-(l2T])-(l22ll does admit a solution, we can obtain this solution denoted by d'^{Qf,Qf), Pi{Qf,Qf), 
P2{QijQf)' ^'^{QtjQf) thanks to a simple procedure inspired by the single-cell procedure proposed 
in Part I [1] for finite number of users: 

• Solve the system ([T9l)- (l20l )-(l2TI)-(l22l ) formed by replacing the equality in equation (l22l) of sys- 
tem (O-dlOll-dlB-dlZl) by the following inequality 

J g-(x, pt, Q\,e)dX'{x) KQi. (EH) 

The existence and the uniqueness of the solution to this new system for an arbitrary {Qf, Qf) € 
can be proved by extending, to the case of infinite number of users. Proposition 1 which was provided 
in III for the case of finite number of users. 
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• If the resulting power 3"(x, Pf, Qi,^'^)d\'^{x) transmitted in the interference band is equal 
to Qi, then the resulting value of d^(Qi^,Qf) coincides with the unique solution to system (fT9l )- 
(|20l)-(l2T])-(l22l) . Once again, this claim can be proved by extending Proposition 1 of HI to the case 
of infinite number of users. 

• If the power 3"(x, P^, Q'l,^^)dX^{x) is less than Q1, then d^{Qf, Qf) is clearly not a solution 
to system (fT9l)-(l20l)-(l2TI)- (l22l) , as equality (l22l ) does not hold. In this case, it can be easily shown 
that system ([T9l)-(l20l)-(l2T]l-(l22ll has no solution. The point {Qf,Qf) is thus eUminated. 

ii) Compute the total power 

QT{Qt,Q?)= E J2^k,lPk,l+li,2Pk,2 
c=A,B k 

that would be transmitted if the values of and Qf introduced by Theorem [T] were respectively equal 
to Qf and Qf. 

Hi) The final value of d^, d^ is given by d^{Qf, Qf), d^{Qf, Qf), the value associated with {Qf, Qf) 
the argument of the minimum power transmitted by the network: 

(Q^^,Qf) = arg min QT{Qt,Qf) ■ 

iv) Finally, we choose 

^^bopt = ^"'i ^Sbopt = d^- 
Note that the same procedure provides as a byproduct the limit Qt of the total transmit power as 

QT = QT(g^Qf). 
Comments 

It is clear from our previous discussion that the above procedure for computing (d^, d^) can be done 
in advance prior to resource allocation. This is essentially due to the fact that the asymptotically optimal 
pair of pivot-distances {d^, d^) does not depend on the particular cell configuration, but on an asymptotic 
or "average" state of the network. The procedure can be run for instance before base stations are brought 
into operation. It can also be done once in a while as the asymptotic distribution of the users and the 
average rate requirement f can be subject to changes: but these changes occur after long periods of time. 
Therefore, the number of operations needed for the computation of {d^,d^) is not a major concern 
because it does not affect the computational complexity of resource allocation. 
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D. Asymptotic Performance of the Simplified Algorithm 

(K) 

Denote by Ql^f^^^ the total power transmitted when our simpUfied allocation algorithm is applied. 
Recall that designates the total power transmitted by the network when the optimal resource 

allocation associated with the joint resource allocation problem (Problem 2 of UJ) is used. 

Proposition 2. The following equality holds: 

Proposition |2] can be proved using the same arguments as the ones used in Appendix |Al The detailed 
proof is omitted. The above Proposition states that the proposed suboptimal algorithm tends to be optimal 
w.r.t. the joint resource allocation problem, as the number of users increases. Therefore, our algorithm 
is at the same time much simpler than the initial optimal resource allocation algorithm of [1], and has 
similar performance at least for a sufficient number of users in each cell. Section IVIII will furthermore 
indicate that even for a moderate number of users, our suboptimal algorithm is actually nearly optimal. 



V. On the Convergence of the Simplified Allocation Algorithm 

As stated before, the simplified algorithm performs the resource allocation in each Cell c independently 
for the protected X'p and the non protected X'j users, which are separated by the predefined pivot-distance 
^subopf Resource allocation for the non protected users is done by the iterative and distributed ping-pong 
algorithm described in Section |llll It was stated in Section |lll] that the convergence of the ping-pong 
algorithm is ensured by the feasibility of the the problem of resource allocation for the non protected 
users {Xf,Xf} (Problem[I]). If Problem [U is feasible, the ping-pong algorithm converges. If Problem [T] 
is infeasible, the ping-pong algorithm diverges. It was also stated in Section |lll] that Problem [1] may not 
be feasible if arbitrary values of the pivot-distances d^bopt ^^'^ '^Sbopt used. Fortunately, feasibility 
of the latter problem will not be an issue if the value of d^bopt ^^'^ ^smbopt relevantly chosen as 
described by the procedure introduced in Section |TV-D[ Indeed, it can be shown in this case that at least 
for large K, the set Xj will contain the users who would anyway have been restricted to the interference 
band J if the optimal resource allocation of Part I HI was used. More precisely, it can be shown that 
there exists a value Kq of K beyond which Problem [J is always feasible. The proof of this statement is 
provided in 1131. It is based on sensitivity analysis of perturbed optimization problems |fT4l . It is worth 
mentioning that in our simulations. Problem [T] was feasible in almost all the settings of the system, even 
for a moderate number of users per cell as small as 25. 
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VI. Selection of the Best Reuse Factor 

The selection of a relevant value a allowing to optimize the network performance is of crucial 
importance as far as cellular network design is concerned. The definition of an optimal reuse factor 
requires however some care. The first intuition would consist in searching for the value of a which 
minimizes the total power = (a) transmitted by the network, for a finite number of users K. 
However, q')^^ {a) depends on the particular target rates and the particular positions of users. In practice, 
the reuse factor should be fixed prior to the resource allocation process and its value should be independent 
of the particular cells configurations. A solution adopted by several works in the literature consists in 
performing system level simulations and choosing the corresponding value of a that results in the best 
average performance. In this context, we cite HI, 191 and ifTOll without being exclusive. In this paper, we 
are interested in providing analytical methods that permit to choose a relevant value of the reuse factor. 
This is why we propose to select the value Oopt of the reuse factor as 

Oopt = arg min lim Q^r^^ (a) . 

Recall that the limiting power = lim/^^oo Q^^^ is given by equation (fTSl) . In practice, we propose to 
compute the value of Qt = Qt{o:) for several values of a on a grid in the interval [0, 1]. For each value 
of a on the grid, Qxict) can be obtained using the procedure presented in subsection IIV-CI Note also 
that complexity issues are of few importance, as the optimization is done prior to the resource allocation 
process. It does not affect the complexity of the global resource allocation procedure. We shall see in 
Section IVIII that significant gains are obtained when using the optimized value of the reuse factor instead 
of an arbitrary value. 

VII. Simulations 

We first begin by presenting the technical parameters of the system model. In our simulations, we 
considered a Free Space Loss model (FSL) characterized by a path loss exponent s = 2 as well as the 
so-called Okumura-Hata (0-H) model for open areas |[T5l with a path loss exponent s = 3. The carrier 
frequency is /o = 2AGHz. At this frequency, path loss in dB is given by pdsix) = 20 logiQ(x) + 100.04 
in the case where s = 2, where x is the distance in kilometers between the BS and the user. In the case 
s = 3, pdB{x) = 301ogxo(3;) + 97.52. The signal bandwidth B is equal to 5 MHz and the thermal noise 
power spectral density is equal to Nq = —170 dBm/Hz. Each cell has a radius D = 500m. 
Asymptotically optimal pivot-distance and frequency reuse factor: We first apply the results of of 
Sections |IV] and |Vl] in order to obtain the values of the asymptotically optimal pivot-distances d^, 
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and the asymptotically optimal reuse factor Oopt- These values are necessary for the implementation of 
the simplified allocation algorithm proposed in Section Unl Each of the two cells is assumed to have in 
the asymptotic regime the same uniform distribution of users: = = A where dX{x) = dx/D. 
The average rate requirement in each cell is assumed to be the same, too: f"^ = = f, where is 
defined in Subsection IIV-A2I as the average data rate in Cell c measured in bits/sec/Hz. In this case, the 
optimal pivot-distance is the same in each cell i.e., = . Define d^pi = d^ = d^ . The value of 
dopt and aopt was obtained using the method depicted by Subsection IIV-CI and Section |Vl] respectively. 
Denote by rt the total data rate of all the users of a sector measured in bits/sec (r^ = f * B). Figure |2] 
and Figure [3] plot respectively Uopt and the normalized pivot-distance dopt/D as functions of the total 
rate rt for two values of the path loss exponent: s = 2 and s = 3. Note from Figure [2] that aopt and dopt 
are both decreasing functions of rt. This result is expected, given that higher values of rt will lead to 
higher transmit powers and consequently to higher levels of interference. More users will need thus to 
be "protected" from the higher interference. For that purpose, the pivot-position must be closer to the 
base station and a larger part of the available bandwidth must be reserved for the protected bands Ta 
and Tb- Note also that, in the case s = 3, "less protection" is needed than in the case where s = 2. In 
other words, dopt(s = 3) > dopt{s = 2) and aopt(s = 3) > aoptls = 2). This observation can be explained 
by the fact that, when the path loss exponent is higher, the interference produced by the adjacent base 
station will undergo more fading than in the case when the path loss exponent is lower. 
Simplified resource allocation: In Section |llll we proposed a suboptimal allocation algorithm charac- 
terized by its reduced computational complexity compared to the optimal allocation algorithm depicted 
in HI- This algorithm assumes fixed pivot-distances d^bopf '^Sbopt- Here, we study the performance of 
this algorithm when d^bopt ^^'^ '^Sbopt chosen according to the procedure provided in Section IIV-CI 
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i.e., c^^bopt = ^opt and d^^opt = ^opt> where dopt is the asymptotically optimal pivot-distance defined 
earlier in this section. In order to study the performance of this algorithm, we need to compare, for a 
large number of system settings, Qsublpt ^^^^^ transmit power that must be spent when applying the 
simplified algorithm, with the total transmit power that must be spent when the optimal resource 

allocation scheme of Part I 1 1 ] is applied. The results must then be averaged in order to obtain performance 
measurements that are independent of the particular system setting. We consider therefore that users in 
each cell are randomly distributed and that the distance separating each user from the base station is 
a random variable with a uniform distribution on the interval [0,D]. On the other hand, we assume 
without restriction that all users have the same target rate, and that the number of users is the same for 
the two cells = . Define x as the vector containing the positions of all the users in the system 
i.e., X = (xi, X2, ■ ■ ■ ,xk'')c=a.b- Recall that V/c, is a random variable with a uniform distribution 
on [0,D]. For each realization of x, define (5^''(x, a) as the total transmit power that results from 
applying the optimal joint resource allocation scheme of Part I with the value of the reuse factor fixed to 
a. Define (^^^(x) = mine (^^^(x, a). In the same way, denote by QsL^opt(^) transmit power 

that results from applying the simplified resource allocation scheme of Section |lll] with the value of the 
reuse factor fixed to Oopt defined in Section |Vll For each realization of the random vector x, the values 
of Q^^(x) and Q^^^^lpii^) were calculated and then averaged to obtain Ex[(5^^(x)] and ]Ex[(5si^opt(-^)] 
respectively. We plot in Figure |4] the values of Ex[Q^^^(x)] and ]Ex[(5si!bopt(^)] ^ range of values 



s=2, D=500m 



^asymptotic 
-optimal (K*=k'^=50) 
-simplified {K*=k'^=50) 
simplified (K*^K°=25) 





Figure 4. Optimal and suboptimal transmit power vs. sum rate Figure 5. Transmit power vs. the pivot-distance d for the 

simplified allocation scheme (rt = 10Mbps, K" — 50) 



of the sum rate rt measured in bits/sec (rt = Ylk=i ^kB) in two cases: K'^ = 25 and K"^ = 50. The 
error bars in the figure represents the variance of the random variable (5si(^opt(x) in the case K"^ = 50. In 
the same figure, the corresponding values of the asymptotic transmit power Qt defined by Theorem [T] 
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are also plotted. This figure shows that, even for a reasonable number of users equal to 25 in each cell, 
the transmit power needed when we apply the suboptimal algorithm is very close to the power needed 
when we apply the optimal resource allocation scheme. The gap between the two powers is of course 
even smaller for K'^ = 50. This result validates Proposition [2] which states that our proposed suboptimal 
resource allocation scheme is asymptotically optimal. Figure [5] is dedicated to illustrate the sensitivity 
of the simplified allocation scheme with respect to the pivot-distance dsubopt in the case K'^ = 50. For 
that sake, the figure plots the total transmit power resulting from applying the simplified scheme as a 
function of dsubopt- The minimum in the figure corresponds to the asymptotically optimal pivot distance 
f^subopt = c^opt- We note that using values different from dopt increases the suboptimality of the simplified 
scheme. Let us go back to Figure ID The latter figure shows that over the range of the considered values 
of the total data rate rt, the total transmit power Ex[Q^^(x)] for K'^ = 50 is practically equal to the 
asymptotic power Qt- This result suggests that, for a number of users equal to 50 in each cell, the 
system is already in its asymptotic regime. In order to validate the latter affirmation, one still needs to 
investigate the value of the mean square error {Q^^ — QtY well. This is done by Figure |6] which 
plots ^''(Qt c^^~^^^ ' mean square error normalized by Qt- 

s=2, D=500, r =5 Mbps 
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Figure 6. (S^ number or users per cell 



VIII. Conclusions 

In this pair of papers, the resource allocation problem for sectorized downlink OFDMA systems has 
been studied in the context of a partial reuse factor a G [0, 1]. In the first part of this work, the general 
solution to the (nonconvex) optimization problem has been provided. It has been proved that the solution 
admits a simple form and that the initial tedious problem reduces to the identification of a limited number 
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of parameters. As a noticeable property, it has been proved that the optimal resource allocation policy 
is "binary": there exists a pivot-distance to the BS such that users who are farther than this distance 
should only modulate protected subcarriers, while closest users should only modulate reused subcarriers. 
A resource allocation algorithm has been also proposed. 

In the second part, we proposed a suboptimal resource allocation algorithm which avoids the costly 
search for parameters such as the optimal pivot-distance. In the proposed procedure, the optimal pivot- 
distance is simply replaced by a fixed value. In order to provide a method to relevantly select this 
fixed pivot-distance, the asymptotic behavior of the optimal resource allocation has been studied as the 
number of users tends to infinity. In the case where the fixed pivot-distance associated with the simplified 
algorithm is chosen to be equal to the asymptotically optimal pivot-distance, it has been shown that our 
simplified resource allocation algorithm is asymptotically equivalent to the optimal one as the number 
of users increases. Simulations proved the relevancy of our algorithm even for a small number of users. 
Using the results of the asymptotic study, the optimal value of the reuse factor has been characterized. It 
is defined as the value of a which minimizes the asymptotic value of the minimum transmit power. Our 
simulations proved that substantial improvements in terms of spectral efficiency can be expected when 
using the relevant value of the reuse factor. 

Appendix A 
Proof of Theorem [T] 

Theorem [T] characterizes the asymptotic behaviour of the minimal transmit power resulting from 
applying the optimal resource allocation when the number of users K tends to infinity. It is thus useful 
at this point to recall the theorem given in the first part of this work which characterizes the optimal 



allocation for finite values of K. Define the function F{x) = E 



For each cell c = A, B and 



for each I = 1 . . . K*^, define by and 6^ the unique positive numbers such that Y^^^^i c{gl\a'=) ~ ^ 
and Y.k=i+i cCg^.Un = with a[j = = by convention. 

Theorem 2 (111). 

(A) Any global solution to the joint resource allocation problem has the following form. For each Cell c, 
there exists an integer € {1, • • • ,K^}, and there exist four positive numbers (31, Qi such that 

1) For each k < L^, 



pc _ ^ mc\-^ f-i I 9k,i{Qi) 
^k,i - 9k,i[Qi) J I i^^c ^'i 

Rk 

7m 



E 



log[l + gk,i{Ql)PkiZ 



(23) 

7^,2 = 
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2) For each k > L^, 



= 



Rk 



(24) 



E 



log 1 + 5^,2^%^ 



3) For k = L'' 



k-l 

«=1 



1 + e 



Pk,2 = 9k,2ni9k,2P'2) 

^ i - a c 
7fc,2 = ^ 2^ 7/,2- 

«=fc+l 



(25) 



(B) For each a = A,B, the system §'^(Q^, Qi) formed by the following four equations is satisfied. 
L = mm <! / = 1 . . . K / F ( ^ , a; ) < gi^2F {gi,2bi) 



1 + e 



7iMC 



+7£^2C(<7L,2/32')=^^L= 



1+e 



(26) 
(27) 
(28) 

(29) 



where the values of ^ and P^^ in i[29\l are the functions of (/3i,/32)C'^) defined by equation f l23l ). 
(C) Furthermore, for each a = A,B and for any arbitrary values and Qf, the system of equations 
§^{Qi,Qf) admits at most one solution (L^, /Jg, '^'^). 



In subsection IIV-A2I we obtained that for each cell c = A,B, 



AT 



r P^/^^ , 0, 0)diy^^^^^ (r, x) + o/<(l) 



(30) 
(31) 



where A^'^^^ = [0,p] x [e,d^'(^)] and Ag'^^^ = [0,p] x [(i^'(^), D] and where ci'='(^) is the pivot- 
distance i.e., the position of user L^'^^\ Our aim is to prove that q'^^ = J2c Qi^^^ + Q2^^^ converges 
as if ^ oo, and to characterize the limit. For each cell c G {A,B}, sequence d^'^^^ is bounded by 
definition {d'^'^^^ £ [0,D]). Consider a subsequence 0/< such that {d'^'^'^'^\d^'^^'^^) converges to a 
certain limit, say {d^,d^). We prove that in this case, all quantities Q^'^'^^^ Q2^''"'\ I3l'^^''\ ()''^^''"'\ 
^c,(0k) converge to some values Q\, Q2, P2' '^^ which we shall characterize. Focus for instance on 
sequence /Jg'^*^^^- Recalling that 72e,(A) 2 tends to zero as K ^ 00 (72e,(K) 2 = ok{^)) and replacing 
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each 7^ ^ with expression (1231 ) 7^ 2 = ^k/C{g'^ 2/^2^^"*)' we obtain immediately 



1 — a 



where we defined 



(32) 



(33) 



for each x, /3, Q, ^. In the asymptotic regime, we obtain the following lemma. 



Lemma 2. As K 00, sequence P'. 



t 



A') 



converges to the unique solution to the following equation: 



D 



rg(x,/32',0,0)di/^(r,x) 



1 — a 



(34) 



Jo Jd" ^ 

Proof: Existence and uniqueness of the solution to (l34l) is straightforward since function /? 1— > 
S(x, /3, Q, ^) is strictly decreasing from 00 to on R_(_. We remark that sequence (i'^'^'^^^ is bounded i.e., 
^c,((^jf) < ^ fQj- ^ certain constant k. In order to prove this claim, assume that there exists a subsequence 
fj'^'^'-^'^^ which converges to infinity. This hypothesis implies that the subsequence given by the Ihs of (|32l ) 
for K of the form K = C,{K') converges to zero as K' 00. This is in contradiction with (|32l ) which 
states that the latter sequence converges to Using similar arguments, it can be shown that ^2^'^'^^ 
is lower bounded by a certain e' > i.e., e' < /Jj'*^'^^'' < k. Denote by P2 any accumulation point of 
^c,((/)k) ^jjJ define ^2^^'^'^ a subsequence of ^2^'^'^'^ (^•^•» (^k coincides with (p(^(^K) for a certain function 
which converges to /32. We prove that P2 is given by (l34l l. Define G{r, x, y) = r^{x, y, 0, 0). We show 
that the difference 

rp N , rp I'D 



A 



10 Jd^-^OK 

tends to zero as ^ 00. By the triangular inequality, 

fP f-D 



A 



K 



< 



Jd' 



P fD 



Gir,x,f32'^'^''^)diy^ir,x) 



+ 



G(r,a;,/32'^^"Vi^''^^"^(^,a;) - / / G(r, x, /Jg'^^" ^^^''^^'H^, ^) 

Jd-'-'.OK) Jo Jd' 

P I'D rp i-D 

/ G{r,x,(52)dv^^'^^'<\r,x) - / / 

Jd" Jo Jd 

P f-D 



Respectively denote by Ak,i, Ak,2, Ak,3 the first, second and third terms of the above equation. We 
first study Ak,i- Clearly, function G{u,x,f3) is bounded on [0, p] x [e,D] x [e',K]. Denote by ^ an 
upper bound. Then, Ak,i < £,1^^^" (Ik), where Ik = [0,p] x d^] (or Ik = [0,p] x [d^ if 

^ (ic,ieK)y Recall that converges to (i^ by definition, so that z^^(/i^) = (''([O, p])A^([d^'(^^), d^]) 
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converges to zero as long as measure A'^ has no mass point at dF. Since z/'^'(^^') converges weakly to 
it is straightforward to show that u^'^^'^^Ik), and thus Ak,i, tend to zero. Now focus on Ak,2- The first 
term J J G{r^x, I32)dv'^''^^'^\r,x) which composes 2 converges to J J G{r,x, l32)di''^{r,x) by the 
weak convergence of u'^'^^^^ to u'^. The second term f f G{r,x,(32^^'^^)dv'^{r,x) converges to the same 
limit by Lebesgue's dominated convergence Theorem. Thus, Ak,2 tends to zero. In order to prove that 



Ak,3 tends to zero, we remark that sup 



9G(x,r,/3) 



< oo, where the supremum is taken w.r.t. {x,r,j3) € 



[0, p] X [e, D\ X [e', k\. Denote by C the latter supremum. We easily obtain G(r, x, fi'^^^'^'') — G{r, x, (32) 



P'/'''^ - P2 



c 

Ak,3 < C 



, so that Ak3<C 



P2 ~ 



< 



u^'^dO, p] X [d*^, D]). Since z^^^ is a probability measure. 



P2 



Thus A;^ 3 tends to zero as K tends to infinity. Putting all pieces together, Ak 



tends to zero. Using (l32l ). | Jq J^o G{r,x,l32 )dv^{r,x) converges to By continuity arguments. 



(32 = limj^ p^^^'^"^ satisfies dH. Thus p^^^'^"^ is a bounded sequence such that any accumulation point 
is equal to /Ja defined by dSI- Thus limj^ (3"^^'''"^ = p2- ■ 

Using Lemma ID we may now characterize the limit of (|3T] ) as ^ oo. Using the fact that 
limx ^2^'^'^^ = ^2 and lim/^ (jpA'f'K) — ^\Q^g with some technical arguments which are similar to the 
ones used in the proof of Lemma |2l we obtain 

j^c fp f-D 



B 



r J(x, 0, 0)dz/^'('^^')(r, x) + ok{1) 



(35) 



where ^2 is the unique solution to (l34l) . As v'^^'^^'^) converges weakly to v'^, Q'^^'^'^^ converges to 



Q'2 



t fp 

'0 



2 .In .Id': 



D 



ry(x,/32,0,0)di/^(r,x) . 



(36) 



The same approach can be used to analyze the behavior of sequences Q^'*^*^^^ and P'^'^'^'^^ for each 
c = A,B. After similar derivations, we obtain the following result. As ^ oo, sequence {pf''^'^^\ 
qMM^ ^AiM^ qB,{^k)^ converges to the unique solution {P^,Qt^'^,P^,Qf,^^) 

to the following system of six equations: 



t 

gi{d'',Ql 



r3^{x,P'l,Ql,e) du%T,x) 



r%x,Pl,Q\,ey%r,x) 



a 



g2{d'')F{g2{d-W2 



A,B, 



(37) 



1 + e V 1 + 

where P2 d'^ are the limits of ^2^'^'^'^ and respectively. We discuss now the existence and 

the uniqueness of the solution to the above system of equation. For that sake, recall the definition of 
functions 3" and S given by ([121) and (l33]l respectively. Note that 3^{x, /3, Q, C) = '3^{x, Q, 0), and that 
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S(x,/3, Q,^) = 9(2;, Q, 0). Define 0^ = for c G {^4,5}. By applying this new notation, The 
first two equations of system (|37] ) give place to the following system of four equations: 



rg(x,/ff,gf,0)i/^(r,2;) =a 



(38) 



The existence and the uniqueness of the solution Qx)c=a.b to the system (l38l ) was thoroughly studied 
in Q. Applying the results of |2| in our context, we conclude that (yrfe) Qi)c=yl,B = (/^fiQi) 



is unique. We turn now back to the third equation of system (|37] ) to get the following equality i^'^ = 

g^(^ii'=)F{g2{d'=)i3'') ^^'^ latter equation proves the uniqueness of for c = vl, B. The uniqueness 

of /3f follows directly from the same equation. 

So far, we have proved the uniqueness of the solution to the system ( [371 ) of equation. As for the 
convergence of sequences (/jf '^'^'^'^ Q^'^^"\ ^^.('^^), f3^'^'^''\ Q^'^'^''\ ^^'('^^)) to this unique solution, 
its proof is omitted here due to the lack of space, but follows the same ideas as the proof of convergence 
of iP^'^'^'^K and '("^"^ Qf •('^"b provided above. 

So far, me managed to prove that for any convergent subsequence (^d'^'^'^'^\d^'^'^^^) — > {d'^,d^), the 
set of parameters {Q''{^''"'\ Q^/''"'\ (i'q^'^'''^ , C''('^'^))c=a,b converges to some values Q\, (31, 

P2, which are completely characterized by the system of equations (l34l ). (l36l ) and (l37l ). as functions 
{d^,d^). Using decomposition = Q'^ y. A^, the system formed by equations ([34]), ( [36] ) and ( [37] ) is 
equivalent to the system ([T9l)-([20l)-([2TI)-([22l) provided in Theorem [T] At this point, we thus proved that at 
least for some subsequences (f)K defined as above, the subsequence Q^^*^^^ converges to a limit which 
has the form given by Theorem [T] The remaining task is to prove that is a convergent sequence. 

First, note that is a bounded sequence. Indeed, is defined as the minimum power that can 
be transmitted by the network to satisfy the rate requirements. By definition, Qj, is thus less than the 
power obtained when using the naive solution which consists in forcing each base station to transmit 
only in the protected band (7^ ^ is forced to zero for each user k of each cell c). Now it can easily be 
shown that when K ^ 00, the power associated with this naive solution converges to a constant. As a 
consequence, one can determine an upper-bound on which does not depend on K. 

Second, assume for instance that Qt and Q'j. are two accumulation points of sequence Q^^. By 
contradiction, assume that Qt < Q'rp. Extract for instance a certain subsequence of which converges 
to Qt- Inside this subsequence, one can further extract a subsequence, say 6k, such that 



qJ?'^) — > Qt, d^, yc = A,B 
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where d,-^ and are some constants both (just use the fact that d^'^^^ is bounded for each c). Clearly, 
Qt can be written as in (fTSl) . where parameters Pf, /Jf, d^, Qi, '^'^ satisfy the system of equations ([T9l )- 
(|20l)-(l2T])-(l22l) . We now consider the following suboptimal resource allocation policy for finite numbers 
of users and . In each cell c G {A, B}, users whose distance Xk to their BS is less than d'^ are 
forced to modulate in the interference band J only, while users k which are farther than d'^ are forced to 
modulate in the protected band Tc only. In other words, for each user k in cell c, we impose 

[C] <^ '''^ '''^ yc = A,B. (39) 

Particular values of the (nonzero) resource allocation parameters "y^i^^ki obtained by minimizing 
the classical joint multicell resource allocation problem (Problem 2 in |1]), only including the additional 
constraint [C]. As a new constraint has been added, it is clear that the total power transmitted by 
the network, say Qj, ^' , is always larger than the total power Qj, achieved by the optimal resource 
allocation, for any K. On the other hand, using the same asymptotic tools as previously, it can be shown 
after some algebra that 

limQ^^^'^'* = limQ^^^'^ = Qt . 

-^l K ^ 

In other words, this suboptimal solution performs as good as the optimal one when K has the form 
K = 9k' for some K' . Although we omit the proof, this observation is rather intuitive. Indeed for such 
K = 6k', the optimal values of the pivot-distances converge to the arbitrary ones d^, d^ . Even more 
importantly, it can be shown that the total power spent when using the suboptimal procedure 



converges as K ^ oo. Therefore, 

limQ^^^'* = Qr 

Now consider a subsequence ij^K such that YmiK Q^^^ = Q'j- > Qt, and compare our suboptimal 
allocation policy to the optimal one for the i^'s of the form K = iI^k'- As \miK Q^t'^'' ^ ImiK Q^'^^'*, 
there exist a certain e > and there exists a certain Kq such that for any K > Kq, 

The above inequality contradicts the fact that Q^'^'^ is the global solution to the joint multicell resource 
allocation problem (Problem 2 in |[T1). Therefore, Q'j, necessarily coincides with Qt- This proves that 
Q^^ converges to Qt- To complete the proof of Theorem \T\ one still needs to prove that for any 
fixed value of {Qf,Qf) G M+, the system formed by equations (fT9l)-(l20l)-(l2TI)-(l22l) admits at most one 
solution. The main ideas of this proof were evoked in the proof of Proposition 1 of [1]. However, the 
complete proof is omitted due to lack of space. 
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